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ABSTRACT 



It is shown that the height scale for accretion disks is a constant whenever 
hydrostatic equilibrium and sub-sonic turbulence regime hold in the disk. In 
order to have a variable height scale, processes that do contribute with an extra 
term to the continuity equation are needed. This makes the viscosity parameter 
much greater in the outer region and much smaller in the inner region. Under 
these circumstances, turbulence is a presumable source of viscosity in the disk. 

Subject headings: accretion disks - hydrodynamics, turbulence, dwarf novae, 
radiative transfer 



1. Introduction 



In the last 35 years, accretion disks have become one of the most intense research ar- 
eas in theoretical astrophysics. The reason for that cannot be explained just by invoking 
their ubiquity in a variety of dif ferent astrophysical environ ment. At that point, one should 
recognize the role played by the IShakura fc Sunyaev! (119731 ) a standard model, which made 
appeal to a very intuitive, simple and neat physics. This model has, as its main assumptions, 
geometrical thinness, large optical thickness, hydrostatic equilibrium in z, perpendicular to 
the plane of the disk, and the viscosity parametrization in terms of an unknown a parameter. 
It was meant to apply under very specific conditions, such as those occurring in the outer 
parts of accretion disks in Bynary Systems, around very You ng Stellar Objects, and in flows 
associated with the central engine for Active Galactic Nuclei (jPapaloizou &: Linlll995l ). Soon, 
people realized that the approach used to treat viscosity was general and could be used under 
quite different conditions. Somehow this brought some ease to the long standing question 



-2- 



of viscosity in accretion disks, encouraging people to tackle related problems, even without 
knowing what the source of viscosity really is. So successful and profuse the applications were 
that, today, one speaks of, at least, four models of accretion disks: the a-standard model, suc- 
cessfully applied to Cataclysmic Variables, Transient X-Ray Sources, acc retion disks around 
Active Galactic Nuclei, a c cretion disks in Young Stellar Obj e cts (YSO) (jSmaklll999l ; lLasota 
200ll : iMenou et all I2OO0I : king et all I2OO7I : iBegelmanl Il985l : Eir] Il989h : the ADAF model, 
where the acronym stands for Advection Dominated Accretion Flows, used to e xplain X and 
7-ray emission from underluminous X-Ray Binaries and Active Ga lactic Nuclei (IPerna et al. 



20001 ; iNaravan fc Yilll995l : Becker fc Lell2003l : iNaravan et all 120021 ); the CDAF model, Con- 



vection Dominated Accretion Flows, applied to underfed and radiat i vely inneficient hard 



X-Ray Binaries and Active Galactic Nuclei (lAbramowicz et al. 



Ieumenshchev 
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Ieumenshchev et al. 


2000; 


Ball et al 


2001 


); the 


Shapiro et al. 



2003; 



Narayan et al 



2002 



2002 



Naravan et al. 



2000 



Quataert &; Gruzinov 



(119761 ) two-temperature accretion disk model ap- 
plied to Cygnus X-l. Despite a remarkable difference among these models and systems to 
which they apply, the a-parametrization works quite w ell. One of the rea sons for that would 
be the weak dependence of the disk properties on a ( jKing et al.l 120071 ). Besides, it seems 
that the spread on the values of a to fit all these systems i s not large, i.e., fairly similar 
values of a result in reasonable agreement with observations (jKing et al.l 120071 ). There are, 
however, some points that require more detailed treatment, and these are related to the disk 
scale height. By this we mean not only its value but, above all, how it behaves along the 
disk . These questions are of fundamental importance when one is concerned with charac- 
teristic time scale lengths, turbulence regimes, and criteria to choose among different energy 
transport models. To make this point more clear, we shall focus on the a-standard model, 
but our criticism applies to them all. The assumed constancy of the viscosity parameter 

g 

is very decisive on the a-standard model, leading to a disk scale height that goes like r», 
and yields the same behavior along r for both the effective temperature and the temper- 

3 

ature at the mid-plane (z = 0) of the disk, T « r~3. As a matter of fact, it leads to a 
constant optical depth along r. This result, apparently consistent, is made possible only by 
an unsound interpretation of a formal solution to the continuity equation. Using the same 
assumptions of the a-standard model, but the constancy of the viscosity parameter, and 
withdrawing the flaw by taking the correct solution to the continuity equation, the results 
change quite significantly. Now, a will go as r 1125 , with a huge variation along the disk. If 
constancy of a is indeed required, one will have to look for some process that add a term 
in the continuity equation in such a way as to give the correct disk scale height. It is our 
claim in this letter that the questions we have just mentioned cast some doubts about the 
way the disk scale height is obtained and its implications in the determination of a. The 
results we have obtained in this letter do modify previous results in the area, and do modify 
our current knowledge about fundamental issues in accretion disk theory. Besides, they are 
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very important as far as full consistency is required and have not been considered so far. 

2. The solution to the continuity equation and the radial height scale profile 

In this section we are going to show how the radial scale height is related to the mass 
continuity equation. We also show that, in order to have a disk scale height varying with 
radial distance, one needs to look for processes that add an extra term to the continuity 
equation. Since the disk is assumed to have azymuthal symmetry, being under hydrostatic 
equilibrium in z-direction, the continuity equation reads, 



with r, p, V r being, respectively, radial distance, mass density and radial velocity. It should be 
said that the above equation neglects mass transport due to turbulence, which is equivalent 
to assuming sub sonic turbulence regime. Now, we set f(r,z) = r pV r and expand / in 
powers of z to obtain 



where account was taken of reflexion symmetry over z. Inserting f(r , z ) given by eq.(2) into 
eq.(l) gives 








(3) 



Setting z = 0, we must have 



d_ 

dr 



/(r0)=0 



(4) 



or 




(5) 



where Co is constant. For z ^ 0, we recall that 





which, by eq.(l), should give 




Q2n 



d 



(7) 



Then 
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and the most general solution to eq.(l) should be written 

n=0 U V ; 

where x = f , and C = (r pV r ) z=0 - Finally, we may write for the accretion rate, 

oo „ 
n=0 v ' 

and since 

^M = 0, (11) 
or 

this implies 

|-^ = °. ( 12 ) 
or 

because the C2„ are all constants. We, then, must conclude that, under hydrostatic equilib- 
rium together with sub sonic turbulence, the height scale of the disk is not allowed to vary 
along the radial distance r. 



3. The dependence of the viscosity parameter on the disk height scale 

In this section, we are going to highlight some consequences of the conclusions we have 
drawn in the last section. In order to proceed, let us recall some results very familiar from 
the accretion disks theory. Let us start from the hydrostatic equilibrium equation, 

^-p = - P n 2 z, (13) 

a z 

which serves to define the disk height scale £, i.e., 

P = pQ 2 £ 2 , (14) 

with P being the pressure, and Q the Keplerian angular velocity. The angular momentum 
conservation equation may be written 

Mttr = 4n-^-(aPl 2 r) , (15) 

from which we obtain 
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a being the viscosity parameter, and 



S 



;i7) 



takes into account null boundary condition for the torque at r = 7*1, the inner radius of the 
disk. Assuming the disk to be cooled by radiative transport in z direction, we may write, in 
the diffusion approximation, 

18 



F = —c P 

with c, i^, r, P r being, respectively, velocity of light, radiative flux in z direction, optical 
depth, and radiation pressure. Replacing differentials by finite differences, and recalling the 
definition of effective temperature, i.e., 



eff 



47T 



Mtt 2 



eq.(18) may be written as 



Ana 



MQ r — T 



(19) 



(20) 



with a being the Stefan-Boltzmann constant, and T the temperature at the mid plane of the 
disk. In the outer parts of the disk, the opacity is mainly given by the free-free opacity. So, 
using a Rosseland mean opacity, averaged over z, eq.(20) will be rewritten as 



>2 2 



T = 2.62 x 10 Zb MVt z p z T 
Finally, for a gas pressure dominated disk, we obtain from eq.(21), 



A>3 22.5 

a^l.SSxlO- 51 ^^, 
M 34 y M 



(21) 



(22) 



where Mn, M 34 , x, y are, respectively, accretion rate in units of 10 g s~ , mass of the 
central object in units of 10 34 g, the radial distance in units of the inner radius, the disk scale 
height in units of the inner radius. The inner radius n is assumed to be 3R g , R g being the 
gravitational radius. According to iKing et al.l (120071 ). if constraints on a are required, one 
should make res ource to observations of system s subject to temporal behavior, such as dwar f 



nova outbursts (IWarner 112003; Cannizzo 



200 



proto stellar accretion disks (jliartmann et al 



, ou tbursts of X-ray transie nts (ILasotall2001l). 



19981 ). FU Orionis outbursts (Lodato fc Clarke 



2QM), or optical variability in Active Galactic Nuclei (jStarling et al.l 120041 ) If we take, for 
instance, the viscous time-scale, 



t = - 



dr , 



(23) 
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where v is the kinematic viscosity, we obtain for the standard model after little algebra, 
using eq.(22), 



_ = £ = 1.63 x 10" 3 17 



225 



™-0.319 y.0.25 _0.125 



r x M 34 a325 d 

x<i being the disk size in our units. Now, inserting this expression into eq.(22), it yields 



(24) 



a s = 207.53 M 34 L375 Mf, - 375 x*' 575 r 125 



or, to make a comparison with ISmak! (119991 ) 



t = 71.39 M34 1 ' 1 M± 7 ' 3 a a ~°' 8 Xd" 25 



(25) 



(26) 



It should be remarked that, contrary to ISmaki (119991 ) paper, we use the disk scale height 
dependence on T, the temperature at z — 0, not T e //, which we believe to be the correct 
procedure. Essentially, this makes ou r scale heigh t greater by a factor of r 0125 , r being the 
optical depth. Our results differ from ISmak! (119991 ) by minor discrepancies in the exponents, 
but by a numerical factor of about two orders of magnitude greater. Finally, using our 
formulation, we obtain 



i = y = ux 10 - ( M " 



0.1875 



r x 



M: 



(27) 
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and for ceo, the value of a directly related to the viscous time, 

a = 0.1x d " L575 a s , 



(28) 



with a s given by eq.(25). Expression (28) is the value of a at x — 1. The value of a anywhere 
in the disk is 

(29) 



From eq.(25) and eq.(28) we see that the viscosity parameter related to the viscous time 
will be much smaller when we employ the correct solution to the continuity equation. At 
the outer radius of the disk it will be 0.1 x^ 9 ' 675 greater than the value obtained with the 
standard model. 



4. Analysis and conclusions 



A rapid inspection on eq.(22) reveals a very strong dependence of a on the radial 
distance, due to the assumptions of sub sonic turbulence and hydrostatic equilibrium. The 
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a parametrization of the viscosity yields hardly credible results. In a disk of Xd = 100, a 
varies by a factor of 3 x 10 22 as compared to its value at x — 1. Assuming equality between 
the length scale of the turbulence and the height scale of the disk, since the disk is assumed 
to be thin, we have for y = 0.01, sub sonic turbulence only for x < 3.04; and for y = 0.1, 
only for x < 23.56. The results are highly dependent on the extent and thinness of the disk. 
The thicker the disk, the more sub sonic it will be. It should be stressed that the conclusions 
we have drawn are based on the analysis of a solution obtained under conditions that, under 
the usual procedure, would result in the a standard model, which gives a disk scale height 

9 

varying with n, a necessary condition needed to have constant a. However, the constancy of 
£ along r is not an assumption, but it stems as a consequence from a rigorous solution to the 
mass continuity equation. In other words, the assumption of constant a is not compatible 
with the solution to the mass continuity equation. To make things compatible, we should 
have 

~(rpV r ) + L(p,v)=0, (30) 

where L = , ^ ) is an operator applied to p and V. If we insist in hydrostatic equi- 

librium, L describes turbulent mass transport. It is beyond our aim, in this letter, to go 
any further on this matter but, if a disk height scale dependence on r is essential to have 
physically meaningful results, a urgent search for physical processes that do contribute with 
an extra term to the mass transport equation is indeed required. In that respect, it is very 
unlikely to discard turbulent mass transport in the disk as one of the reasons to have the 
disk scale height varying along the radial distance, and, if so, turbulence is a presumable 
source of viscosity in the disk. The points we have raised in this letter deserve attention by 
themselves but in no way exhaust the subject. The value of the disk height scale and the 
way it behaves along r are decisive to know, among other things, what is the relevant energy 
transport mechanism in a given region and how it varies along the disk. It is our intention 
to to consider these issues, in a more detailed way, in a future article. 

The authors acknowledge support from FAPESP . 
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